Interface Landau levels in graphene monolayer-bilayer junction 
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Electronic structure of graphene monolayer-bilayer junction in a magnetic field is studied within 
an effective-mass approximation. The energy spectrum is characterized by interface Landau levels, 
i.e., the locally flat bands appearing near the boundary region, resulting in a series of characteristic 
peaks in the local density of states. Their energies are independent of boundary types such as zigzag 
or armchair. In the atomic scale, the local density of states shows a Kekule pattern due to the valley 
mixing in the armchair boundary, while does not in the zigzag boundary. 



I. INTRODUCTION 

Graphene 1 -"— and its bilayer— ~— are characterized by 
zero-gap band structures supporting different types of 
chiral particles and Landau-level structures. Recently, 
atomically thin graphene samples were experimentally 
fabricated using mechanical exfoliatio n 17 i 18 and epitax- 
ial growth. 19,2 - The characteristic Landau- level structure 
and integer quantum Hall effect 1 "^- 7 - i 9 -^ were observed in 
magnetotransport measurements ! 21 ' 22 In this paper we 
study the electronic structure of a hybrid system com- 
posed of monolayer and bilayer graphenes in magnetic 
fields. 

The band structure of monolayer graphene is charac- 
terized by Dirac-like spectrum in which conduction and 
valence bands with linear dispersion stick at the K and 
K' points located at a Brillouin zone corner^— & which 
are called valleys. Bilayer graphene has a zero-gap struc- 
ture, but with quadratic dispersion unlike monolayer.^— 
In a magnetic field, the level structure of monolayeri^— 
and bilaye r 9 ' 15 ' 23 ! 24 differs in number of degeneracy at 
zero-energy and the quantum Hall plateaus appear at 
different filling factors accordingly] 15 i 21 i 22 

The electronic states of graphene with an edge have 
been studied in theories^ - — In particular, when the 
boundary is along zigzag direction, special states local- 
ized at the edge appear as zero-energy modesi 25 ' 26 Similar 
zero-energy edge states exist also in bilayer graphenei 46 ' 47 
In a magnetic field, electron-like and hole-like Lan- 
dau levels are shifted upwards and downwards near the 
boundary, respectively, forming edge channels away from 
zero energy 15 i 29 i 48 ~ — Recently, the transport through 
quantum structures consisting of monolayer and bilayer 
graphenes was investigated] 51 ' 52 In a previous paper, 
the boundary condition between monolayer and bilayer 
graphenes connected by a monoatomic step was studied, 
and the transmission probability through the junction 
was calculated in the absence of magnetic field^ 

In this paper, we study the energy spectrum and lo- 
cal density of states of the monolayer-bilayer graphene 
junction in magnetic fields. Based on the previous 
study, 53 * we consider a composed system of half-infinite 
graphene monolayer and bilayer connected by a mono- 



atomic step along zigzag or armchair direction. In Sec. 
ILTl we present effective mass description for monolayer 
and bilayer graphenes and introduce formulation to de- 
scribe Landau levels of the junction in Sec. IIIII In Sec. 
IIV1 we numerically calculate the energy spectra for sev- 
eral types of the boundaries as well as the local density 
of states. The conclusion is presented in Sec. [V] 



II. EFFECTIVE MASS HAMILTONIAN 
A. Monolayer graphene 

Graphene is composed of a honeycomb network of car- 
bon atoms, containing a pair of sublattices, denoted by 
A and B. Electronic states in the vicinity of K and K' 
points in the Brillouin zone are well described by enve- 
lope functions (Fj£,Fg) and (Fj£ ,F§ ), respectively, 
in an effective-mass approximation. At the K point, the 
effective Hamiltonian for (Fj£,Fg) isir-^ 
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where v sa 1 x 10 6 m/s is the band velocity ir± — Tr x ±iir y , 
and TV = — z/iV + (e/c)A with vector potential in the 
Landau gauge, A = (0,Bx), giving external magnetic 
field B = V x A. The Hamiltonian at the K' point is 
obtained by exchanging tt± in Eq. (JTJ. 

The wavenumber k y remains a good quantum num- 
ber in the present geometry. The operator 7r± can be 
expressed as 



vir + = ihuiBCL , 



(2) 



where Hub = VZhv/ls with magnetic length Ib = 
\J ch/(eB) and and a are raising and lowering opera 
tors, respectively, defined by 
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Here, the center coordinate of the cyclotron motion is 
defined by 
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k I 2 



The Schrodinger equation then becomes 
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The independent solutions of Eq. (JT)) are given by D v (z) 
and Z?_j / _ 1 (— iz), where D u (z) is Weber's parabolic 
cylinder function defined by 
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with F(a;, 7;z) being Rummer's hypergeometric func- 
tion. The components F£ and F^? are related by Eq. 
d6l) with formula 
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D v (z) and D v (—z) can also be chosen as independent so- 
lutions of Eq. (J7J|, as long as l/r(— z/) is nonzero, i.e., v 
is not or a positive integer. Then, D v (z) and D v (—z) 
exponentially diverge in limits z = — oo and +oo, respec- 
tively, while converge to zero in the opposite side. They 
can never be a bulk eigenfunction, but may appear when 
the system is half-infinite in the x direction. For a non- 
negative integer n, D n (z) and D n (—z) are linearly de- 
pendent, and coincide with usual Landau- level function 
except for a normalization factor as 

D n (z) = (~l) n D n (-z) = 2- n / 2 e- z2 / i H n (z/V2), (12) 

with Hermite polynomial H n {z). The other solution 
D- n -i(iz) then diverges both in limits z = +oo and — oo 
and is excluded. D n (z) at a negative integer n generally 
diverges for z — > — oo. At n = — 1, for example, we have 



F>_i(z) = v^e 2 
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FIG. 1: Plots of u(e), v+(s), and f-(s) with energy e set to 
the vertical axis. Black and white circles represent the Landau 
levels of bulk monolayer and bilayer, respectively. = 3 

is taken for bilayer. 



with error function 



Let us define 



erf (a;) = 



dt. 
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where L and R represent the solutions finite in limits z — ¥ 
—00 and +00, respectively. We will consider a monolayer- 
bilayer junction in which the region x < is monolayer 
and x > is bilayer. The eigen function in monolayer is 
given by 

_ ( iai<f>u- 
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The wavefunction at the K' point can be obtained by 
(Ff,Ff) = (F B <,F«). 



The Landau level energies of bulk monolayer graphene 
are given by the condition that the wavefunction is finite 
in limits x = ±00, i.e., v(e) is non-negative integer n. 
We gelMi 
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The plot of v{e) and bulk Landau-level energies are 
shown in Fig. [T] 



B. Bilayer graphene 

Bilayer graphene is a pair of graphene layers arranged 
in AB (Bernal) stacking and includes A-y and B\ atoms 
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on layer 1 and A 2 and B 2 on layer 2. The layers are 
arranged such that sites B 1 and A 2 are directly below 
and above each other, which are connected by interlayer 
coupling 71 ~ 0.39 eV£l The effective Hamiltonian for 



K 



AV F&, F« 2 , F§ 2 ) is given by^ 



We assume that the system infinitely continues in the 
y direction parallel to the boundary. The wavefunctions 
of monolayer and bilayer regions are required to converge 
in limits x — — 00 and 00, respectively. At given energy 
e, they can be written for the monolayer part (x < 0) 
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The Hamiltonian at the K' point is obtained by exchang- 
ing 7r± in Eq. (fT9|) . 

The eigenfunction of Eq. ([19")) finite in limit x — > +00 
is written as 
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where \i = ± is another degree of freedom, e — 
s/{hu)B), and 71 = 71/ '(TujJb)- The wavefunction at 
the K' point is obtained by (F% , F§[ , Pf 2 ' , Ff 2 ') = 

, i. K 1. K jpK 1. K : 
\ t B2i t A2i t Bn t Al>- 

The Landau levels of bulk bilayer graphene are ob- 
tained by the condition that the wavefunction of Eq. 
(l2Cfl) is finite in limits x — > ±00, i.e., includes only 
of non- negative integer n. Allowed indexes are v+{£) — 
0,1,2,- •• and V-{e) = -1,0,1, For z/_ = -1,0 and 
v + = 0, the wavefunction Eq. (l20l) appears to include 
(j)^ with negative n, but corresponding coefficient such 
as vanishes. At e = 0, there are two energy levels 
for v + = and V- = — l. 9 Figure [T] illustrates v±{e) and 
bulk Landau-level energies of bilayer graphene. 



III. MONOLAYER-BILAYER JUNCTION 

We consider a composite system of monolayer and bi- 
layer graphenes, where the left half (x < 0) is monolayer 
and the right half (x > 0) is AB-stacked bilayer. We 
assume that one layer of the bilayer part, containing A\ 
and B\ sites, seamlessly continues to the monolayer part 
with A and B sites, while the other layer composed of 
A 2 and B 2 sites is sharply cut at the boundary chosen as 
x = 0. In the following, we consider two kinds of zigzag 
boundaries, zigzag-1 (ZZ1) and zigzag-2 (ZZ2) and arm- 
chair boundary (AC) as illustrated in Fig. [5] (a), (b), and 
(c)£L 
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and for the bilayer part (x > 0) 



fI 

* A2 

\F§2j 

b A2 



J 



-iXy/ll 



-iXy/ll 



(25) 



(26) 



with six unknown coefficients A K , A K , B±, and B± to 
be determined by the specific boundary condition. 



A. Zigzag boundary, ZZ1 

The boundary ZZ1 is parallel to the zigzag direction of 
honeycomb lattice, and the front-most line of bilayer part 
is formed by B\ and A 2 sites. As the zigzag boundary 
does not mix the wavefunctions at different valleys K 
and K 1 , the boundary condition is separately expressed 
for each valley. The conditions are^ 



F v M (0,y) = Fl(0,y), 

Fgtfrv) = n(P,y), 

F£ 2 (0,y) = 0, 



(27) 



for v = K and K' . For the K point, the conditions are 
rewritten with use of the wavefunctions (f2"3"|) and (|25[) as 




with 
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where the wavefunctions such as (j>„ represent the values 
at x — 0. For each X, the eigen energies are obtained by 
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(a) Zigzag-1 (ZZ1 




(b) Zigzag-2 (ZZ2) 

■J 




where only the third condition is different from Eq. (|28l) . 
Similarly to ZZ1, we obtain the matrix for the K and K 1 
points 



1 P2<t>? h (33) 



-a 2 (, 



% o'J, 



1 Pttil % oj! 



(34) 







C. Armchair boundary (AC) 

The boundary conditions for the armchair boundary 
AC aviM 



FZ 1 (p,y) = FX(0,v), 



b A2 
F B2 



0, 
0. 



(35) 



where the third and fourth conditions mix the wavefunc- 
tions at the K and K' points. They are rewritten as 
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(36) 



FIG. 2: Monolayer-bilayer graphene junctions with boundary 
types of (a) ZZ1, (b) ZZ2, and (c) AC. 



searching for solutions of det M ZZ1 — 0. The correspond- 
ing equation for K' is 



with 
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D. Interface Landau levels 



(37) 



(38) 



(39) 



B. Zigzag boundary, ZZ2 

The boundary ZZ2 is another zigzag boundary where 
the front-most line of bilayer part is formed by B 2 sites. 
The boundary conditions are^ 



F v A2 (0,y) = 0, 



(32) 



Let us consider a special state of valley v, which satis- 
fies the conditions 



FX 1 (0,y) = Fl(0,y), 

^Bi(o ) y) = ^(o,y), 

i^ 2 (0,y) = 0, 
Fg 2 (0,y) = 0. 



(40) 



Because these include both boundary conditions for ZZ1 
and ZZ2, such a state must be shared by both ZZ1 and 
ZZ2. Those states exist at different series of points (e, X) 
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for v — K and K' , denoted by Pk and Pk>, respectively. 
Further, the wavefunction satisfying Eq. (|4TJ)) at valley 
v also meets conditions (|35|) for the armchair boundary 
when the wave amplitudes of the other valley (opposite 
valley of v) are all zero. As a result, points Pr and Pk* 
are also shared by Landau levels in an armchair bound- 
ary. 

Using some algebra, we can show that at the special 
points Pk and Pk> , the gradient of the Landau level en- 
ergy in X vanishes in any types of boundaries ZZ1, ZZ2, 
and AC. We can show that the second derivative also 
vanishes for ZZ2. The detailed proof is presented in Ap- 
pendix El Accordingly the density of states diverges at 
the identical energies independent of the boundary type. 
Further, at those points, the wavefunctions of monolayer 
part and bilayer part connect smoothly on layer 1, be- 
cause the amplitude on layer 2 locally vanishes and thus 
hardly affects the electron motion on layer 1. As a result, 
the wavefunctions on the monolayer and bilayer sides are 
coupled well, and the amplitude is almost equally dis- 
tributed to both sides. 

As will be demonstrated in numerical results presented 
in the next section, in ZZ1 and AC, e(X) takes a lo- 
cal maximum at each Pk and Pk> in positive energies, 
while there usually exists another point nearby where 
e{X) takes a local minimum, giving divergent density of 
states as well. We will show that, around these points, a 
crossover takes place from a monolayer edge-state mainly 
localized in monolayer, to a bilayer edge-state mainly lo- 
calized in bilayer, when X is varied. It is natural that 
slight shift in X does not change the energy at such 
crossover points, because they are anti-crossing points 
between intersecting energy levels of monolayer and bi- 
layer edgestates. In ZZ2, the energy minima and maxima 
are degenerate corresponding to vanishing second deriva- 
tive, and thus e(X) is even smoother and the divergence 
in the density of states is stronger than ZZ1 and AC. 
These nearly flat-band regions around extrema of e(X) 
can be referred to as the interface Landau levels. 



E. Zero energy levels 

The energy spectrum of a monolayer-bilayer junction 
approaches that of bulk monolayer and bilayer graphenes 
in the limit of X — > +oo and — oo, respectively, be- 
cause the wave function, centered at i = X, mostly 
resides in the bulk region far from the boundary. On 
the other hand, the zero-energy level is special in that it 
is contributed not only by the bulk Landau levels, but 
also by the zero-energy edge states, which are localized 
near the boundary region on the terminated layer of bi- 
layer graphene.— We can analytically obtain the ener- 
gies and wavefunctions of zero energy Landau levels us- 
ing the above formulation, as demonstrated in Appendix 
Efor ZZ2 boundary. Table U summarizes the degeneracy 
of zero energy levels in the limit of A — > ±oo for each 
boundary type, where +1 represents the additional de- 
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(c) AC 
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TABLE I: Number of zero-energy Landau levels per spin in 
the limit of X = ±oo, for (a) zigzag-1, (b) zigzag-2, and (c) 
armchair boundaries. +1 represents extra degeneracy due to 
the zero-energy edge mode. 

generacy due to the edge states. In ZZ1 and ZZ2, the 
edge state appears either of X ~ ±oo depending on val- 
leys, while it is absent in ACj^ 

F. Local density of states 

In monolayer graphene, the amplitude of the wavefunc- 
tions at A and B sites are written in terms of effective- 
mass envelope functions as 8 

Vm(R) = e lKR Ff (R) + e 4 "e tK '- R Ff' (R), 
^ fl (R) = —u}e iv e iK ' R Fg (R) + e iK '- R Ff' (R), 

where rj is the angle between the x axis and zigzag di- 
rection of honeycomb lattice and u) = e 27 ™/ 3 . In bilayer 
graphene, the amplitude can be written as-^ 3 

Vui(R)= e iK - R Ff x (R) + e^e«'- R F£'(R), 
Vm(R)= -c^V^Fj^R) + e iK '- R ^f;(R), 
VM 2 (R)= -^ e "'e lK ' R Ff 2 (R) + e iK ' R Ff 2 ' (R), 
V> B2 (R) = w-V^V^Fj^fR) + e- i,7 e iK '- R Fj 2 '(R). 

(42) 

In a tight-binding model, the local density of state 
(LDOS) on site A located at the position R is defined 
by 

PA {e- R) = J2 ^ " £ (Q) )l^ a) (R)| 2 , (43) 

a 

where and tp^ are the energy and the wavefunction 
of eigenstate a. Similar expressions can be written down 
for other sites B, Al, etc. When the wave amplitudes at 
the K and K' points coexist in a single eigenstate, LDOS 
has a Kekule pattern due to the interference between 
the factors e lKR and e lK R ^ In the present case, this 
is expected to appear in the armchair boundary which 
mixes the K and K' valleys, while absent in ZZ1 or ZZ2, 
where every eigenstate is a single-valley state. We also 
define the spatially averaged LDOS for site A as 

P J(e;R)=J2 £ S(e-e^)\(F^y A (R)f. (44) 

a v=K,K' 

This is an average of the original LDOS in Eq. (|43|) over 
several unit cells in the region smaller than typical length 
scales of the envelope function. 



6 




FIG. 3: Energy spectrum as a function of X in boundary (a) ZZ1, (b) ZZ2, and (c) AC, at magnetic field ?wjb = 7i/3 
(B ~ 10T). Black and white circles represent Pk and Pk> , respectively. Numbers between the levels indicate bulk filling factor 
in limits X = ±00. 



IV. NUMERICAL RESULTS 

Figure [3] shows the energy spectra against X, numer- 
ically calculated for the junctions of ZZ1, ZZ2, and AC 
boundaries at magnetic field of Hojb — 7i/3 {B ~ 10T). 
Landau levels approach those of bulk monolayer and of 
bilayer in the limit X — > 00 and —00, respectively. In 
the boundary region, the valley-degenerate levels split 
and connect to different levels in the opposite side. The 
black and white circles represent the interface Landau 
levels Pk and Pk', respectively, which are independent 
of boundary type and corresponds to local band maxima. 
In accordance with the argument in Sec. IIIIDi we actu- 
ally see that energy levels pass through those points in all 
three cases and the gradient vanishes there. In ZZ1 and 
AC, the band minima are also present near the maxima 
at Pk and P' Kl while in ZZ2 the minima and the max- 
ima merge into inflection points as the second derivative 
vanishes. 

The oscillatory band structures appearing in the 
boundary region can be understood in relation to ter- 
minated monolayer and bilayer graphenes. Let us take 
ZZ1 boundary, and consider a system with infinite on-site 
potential added on an array of B sites near the bound- 
ary, as illustrated as white circles in the top panel of 
Fig. [U The system is then separated into monolayer ter- 
minated with Klein's edge and bilayer terminated with 
zigzag edge. In the effective mass approximation, this is 
equivalent to the boundary condition Fg = for mono- 
layer and Fg X = Fg 2 = for bilayer. Note that in the 
effective-mass approximation, shifting of on-site poten- 



tial position by the order of the lattice constant does not 
make a difference in the result. 

Lower panels of Fig.Ushow the energy spectrum of ZZ1 
(solid lines) and that of the terminated system (dashed 
lines), for each of K and K' . In the terminated sys- 
tem, the independent Landau levels of monolayer and 
bilayer sharply go up as X goes over the boundary i 29 i 48 
Apparently, the spectrum of ZZ1 resembles that of the 
terminated system, with an energy gap opened at ev- 
ery crossing point. The resemblance of the two different 
spectra may be attributed to following reasons. In the 
monolaycr-bilayer junction, when a low-energy electron 
travels from the monolayer to the bilayer, it feels as if B 
sites suddenly disappear at the boundary, because in bi- 
layer, B\ is coupled with A2 to make high-energy states 
away from e — 0. This effect can be roughly modeled 
by on-site potential at B sites at the boundary. For an 
electron coming from bilayer side, on the other hand, 
B\ site which was absent in the low-energy spectrum 
suddenly resumes at the beginning of the monolayer re- 
gion, while A\ just smoothly connects to A. This should 
roughly correspond to some condition for B\ sites, with 
A\ left intact. Energy gaps opening at crossing points 
are due to finite hybridization between monolayer and 
bilayer states. 

As another remark, we observe that energy levels of 
ZZ1 pass through every crossing point of terminated bi- 
layer and monolayer levels. This occurs when an eigen- 
function of ZZ1 happens to have a node on the on-site 
potential sites, because such a state is also an eigenstate 
when on-site potential is present. Therefore the wave- 
function of ZZ1 becomes identical with that of the ter- 
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Zigzag-1 




Zigzag-2 




FIG. 4: (Above) Monolayer-bilayer junction of type ZZ1, sep- 
arated by infinite potential on white circles into independent 
monolayer and bilayer graphenes. (Below) Energy spectrum 
of junction ZZ1 (solid) and separated system (dashed), in 
magnetic field ftwB = 7i/3. Left and right panels show the 
spectra of the K and K' points, respectively, and the black 
and white circles are the interface Landau levels. 



minated system at each crossing point. 

Similar analysis is also available in boundary ZZ2. Fig- 
ure [5] compares the energy spectrum of ZZ2, and that of 
separated system with on-site potential on B sites illus- 
trated in the top panel. The boundary condition be- 
comes Fg = for monolayer and Fg-y = F^ 9 = for 
bilayer. Since the low-energy spectrum of the bilayer is 
dominated by Al and B2, the second condition F V A2 = 
should give a weaker effect compared to F B2 = in ZZ1 
and thus leads to better coupling between the monolayer 
and bilayer region. In Fig.O indeed, the mixing between 
terminated levels looks stronger than in ZZ1, resulting in 
the monotonic dependence rather than non-monotonic 
behavior in ZZ1. 

Figure [5] plots the wavefunctions near an interface Lan- 
dau levels associated with the K' point in ZZ1, where (a) 
to (f) correspond to the points in the energy spectrum in 
Fig. |U The point (b) is exactly at a local maximum Pjc' ■ 
There, the wavefunction of layer 1 smoothly connects at 
the boundary as argued in Sec. 11111)1 while generally not 
in other cases. The point (e) is exactly at the cross- 
ing point of terminated levels. There, the wave function 
indeed has a node at the interface for the B and Bl 



FIG. 5: Plots similar to Fig. |U for ZZ2. 



components and thus can be an eigenstate of the sepa- 
rate monolayer and bilayer. At the local band minimum 
(d), the wave function does not have special features in 
contrast to Pk'- 

In the energy spectrum, the region between (a) and 
(b) and between (e) and (f) have slope close to that of 
the terminated bilayer, while the region between (b) and 
(d) has slope close to that of the terminated monolayer. 
Correspondingly, the wavefunctions of (a) and (f) has 
significant amplitudes in bilayer side, while (c) has great 
amplitude in monolayer side. 

Figure [7] illustrates the typical atomic-scale LDOS of 
Eq. (l4"3")l at interface Landau levels. We here take the 
boundary ZZ2 at the energy of the lowest interface Lan- 
dau level near e = 0.631fiws in Fig. [3^b). The ratio of the 
magnetic length to the lattice constant, Ib/o,, is about 30 
at this magnetic field. The areas of circles in upper and 
lower panels represent the relative amplitude of LDOS at 
each atom on layers 1 and 2, respectively, while open and 
filled circles represent the A and B sublattices. The re- 
sult mainly reflects the wavefunction of interface Landau 
level, since the flat band gives a dominant contribution 
to LDOS. We see that the wave amplitude on layer 1 con- 
nects smoothly at the boundary region, as the amplitude 
of layer 2 is almost absent there. 

For comparison, we show the similar plot of LDOS of 
the armchair boundary at a different energy e = O.SHujb 
in Fig. [S] In accordance with the previous argument, the 
plot clearly exhibits the Kekule pattern unlike in ZZ2. 
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FIG. 7: Local density of states of boundary ZZ2 with fkdB — 
Landau level of K. The areas of circles in upper and lower panels 
and 2, respectively, while open and filled circles represent A and 
field. 



7i/3, at the energy e = 0.63ffius near the lowest interface 
represent the amplitude of LDOS at each atom on the layer f 
B sublattices, respectively. Is jo, is about 30 at this magnetic 



Note that, even in the armchair boundary, the Kekule 
pattern disappears when the energy comes to an interface 
Landau level, because the eigenfunction becomes a single- 
valley state there. 

Figure [5] shows the averaged local density of states 
(LDOS) defined in Eq. flUD, for ZZ2 boundary. The ver- 
tical scale is shared with the corresponding energy spec- 
trum at left. We observe series of peaks corresponding to 
the interface Landau levels of Pk and Pk 1 owing to the 
large LDOS due to the flat band, and its spatial distribu- 
tion is characterized by node pattern of the correspond- 
ing wavefunction. While not shown, the peak patterns 
are quite similar among ZZ1, ZZ2, and AC, since every 
interface Landau level appears at the identical energy 
with the identical effective-mass wavefunction. In ZZ1 
and AC, the band minima appearing near Pk and Pr' 
also contribute to the LDOS divergence and the peak 
structure is a little blurred. 

Near the interface Landau levels, LDOS has a consid- 
erable amplitude in monolayer region, while otherwise it 
is localized mostly in the bilayer region. This is because 
the monolayer and bilayer states are well hybridized near 



the interface Landau levels, while in other regions where 
the band lines are downslope, the states mainly originate 
from bilayer, as argued above. 



V. CONCLUSION 

We have studied electronic structures of monolayer- 
bilayer graphene junctions in magnetic fields. The energy 
spectrum near the boundary region is characterized by 
the interface Landau levels where the band energy is lo- 
cally constant, which arise from hybridization of Landau 
levels of terminated monolayer and bilayer graphenes. 
The energies of interface Landau levels are insensitive to 
the way the second layer is terminated, suggesting that 
they would be robust even in a disordered junction con- 
taining a random atomic configuration at the boundary. 
Interface Landau levels give a characteristic peak pattern 
to LDOS, which may be observed by scanning spectro- 
scopic measurement 
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Here we prove that, at a special point (s,X) where a 
wavefunction satisfies Eq. (|40)) , the derivative of Landau 
level in X vanishes in ZZ1, ZZ2, and AC, and the second 
derivative also vanishes in ZZ2. The condition Eq. PU|) 
is alternatively written as 

det M| Z1 = det M Z ' Z2 = 0. (Al) 

In the following, we will show that Eq. (|A1|) leads to 

— det M Z V = — det M A c = (A2) 

— det Ml Z2 = ^ det Af z u Z2 = 0, (A3) 
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which immediately proves the statements above. 

In the matrices M ZZ1 and M ZZ2 , we can eliminate <j> v '_i 
by replacing them with <fi^' R and 4>^!^[, using the recur- 
sion formula of Weber's function, 

D u+1 - zD v + vD v -i = 0. (A4) 

Eq. (|A1[) at v = K can then be transformed as 

£ Ke 2 ~ + #+i<) = 0, (A5) 

X /i( £ 2 - y fi )4>% + MK+i + = °> ( A6 ) 

At=± 



with /i = — /j,, leading to 



0, 



(A7) 



for both of /i = ±. 

The derivative of the matrix determinant in X can be 
evaluated using Eq. (flO|) . For M ZZ1 , we obtain 



ddetM^ 



dX 



"p)« +1 > 



) 

(A8) 



where ' represents the derivative in A. At the points 
satisfying Eq. (|A1|) . the first term becomes zero because 
of Eq. (|A7|I . The second term is transformed with Eq. 
~T7b as 



^+i ^ 2 



(e 2 - ^)(e 2 



(A9) 



which vanishes since the argument inside the summation 
is anti-symmetric in \i. 



For M ZZ2 , we have 



dX 



_2 fJ.Lj.R 



0^+1+^+1^), 

(A10) 



which similarly vanishes under the condition Eq. (|A7|) 
Eq. (|A10[) is even differentiated as 



E M=± m(£ 2 - ^ 2 - (AH) 



dX 2 



The first term becomes zero again under Eq. (|A7I) . The 
factor )' in the second term can be written in terms 
<f>^_ and (j>^. + i using Eq. (|10p . Then it is shown to vanish 
by similar transformation to Eq. (|A9|) . 

The determinant of Mac can be written in terms of 
those of ZZ1 and ZZ2 as, 

det Mac = det M ZZ1 det M||' 2 - det M ZZ2 det M^ . 

(A12) 



K 
ZZ2 



Under the condition det Mzzi = det M-___ 
(detM^)' = (detM^)' = 0, Eq. (jAl~2l) immediately 
gives (det Mac)' = 0. 



Appendix B: Nearly-zero energy states 

Let us focus on the eigenstates in the vicinity of zero 
energy, taking the case of ZZ2 as an example. We 
will show here that the zero-energy levels in monolayer- 
bilayer junction arc contributed not only by the Lan- 
dau levels in bulk monolayer and bilayer, but also by the 
zero-energy edge states localized to the boundary. For K- 
point, there are two independent wavefunctions exactly 
at zero energy, 



Kl 



A 



1 A2 




(x < 0); 



(x > 0), 



(Bl) 



and 



K2 



3 




: 












°\ 









F K 

r A2 







\F^J 




Wo/ 



(x < 0); 
(x > 0), 



(B2) 



with cj) n = (fi^ = (— l) n 0^ for a nonnegative integer n, 
and the overall normalization factor is omitted. 

In X — > oo, i.e., when the center coordinate goes deep 
inside of the bilayer region, ty K1 and ^> K 2 approach the 
wavefunctions of two zero-energy Landau levels of bulk 
bilayer graphene. In X — > — oo, on the other hand, 
becomes the only zero-energy level of the bulk monolayer, 
while fy K2 does not have any amplitudes in the mono- 
layer side, but mostly concentrated on Bi sites near the 
boundary. ^ K2 at B2 then approximates 



const, x e 



(B3) 



which is independent of magnetic field. This corresponds 
to the zero-energy edge mode in zero magnetic field 
limit 1^ 

For if'-point, we have a single state at zero energy, 



K'l 







1) 


3 






F^\ 




( <h\ 


F& 







1 A2 







\F%2) 




V / 



[x < 0); 



{x > 0). 



(B4) 



When X moves from —00 to +00, the wavefunction ^ K 1 
crosses over from the only zero-energy level in monolayer 
to one of zero-energy levels in bilayer, iv_ = — 1. 
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Besides, for positive large X, there exist another two 
levels near zero energy, which are expressed as a hy- 
bridization of bilayer's Landau level of v + = and a 
zero-energy edge state. The derivation goes as follows. 
By expanding v± in Eq. (|2ip in e, the determinant of 
M222 can be written in a small \e\ as, 

detM|£ a = --<(<^e. + ^- x €.+i) + 0(e 3 ). 
7i 

(B5) 

The energies of the nearly zero-energy states in ques- 
tion are given by the condition tjfi — 0, when Eq. (|B5[) 

vanishes. For a large X, the function <f>^ can be evalu- 
ated by the asymptotic expansion of D v (z) which stands 
for large \z\: 

( Dl 1] (z) (\axgz\ <3tt/4) 
A,(«) ~ < d£\z) + e^D^iz) (B6) 
[ (7r/4 < ± argz < 57r/4) 

!#>(*) = e-VV E(-l) fc ^ " !) 2fc + ^ 



fc=0 



er/* z -»- 1 



k\2 k z 2k 



(B7) 



r(-iz) 

x ^ ( I / + i)(i/ + 2)..-(i/ + 2A) 



fe=0 



fc!2 fc z 2fc 



(B8) 



When z is negative and is small, it approximates 



(B9) 



This leads to an approximate expression </>^ + (x = 0) for 
positive large X, 



b R 



e -wi B y/2_ v ^ 



e (x/i B y/2 

V2X/Ib ' 



(BIO) 



0y becomes zero at v + = (X/ls)e 
the energies of nearly-zero energy mode, 



/7T, giving 



.K'2 



x/h 



■7? 



. e -Wi B ) 3 



(Bll) 



where we used v+(e) ~ (1 + 7i)e for small £. The cor- 
responding wavefunctions in the bilayer part are written 

as 



K'2 



*bulk - 



T bulk — 



K'2 



jK'2 


edge 


^ + ^ 




+ /7i 











1^100- 








(B12) 



(B13) 



(B14) 



V*(l + 7i 2 )(|£|/7i)0-i+, + / 
In X — > 00, the energy 2 becomes exponentially small 

and ^b u i k coincides with the zero-energy Landau level 
of bilayer, v+ = 0. For ^^g 2 ,,, B2 component is nearly 
proportional to and approximates oc e~ kyX near 

x = 0. This is a zero-energy edge state localized near the 
boundary in the bilayer region^ Thus the states 2 
are described as a hybridization of the bulk bilayer Lan- 
dau level and zero-energy edge states. 



1 J. W. McClure, Phys. Rev. 104, 666 (1956). 

2 J. C. Slonczewski and P. R. Weiss, Phys. Rev. 109, 272 
(1958). 

3 D. P. DiVincenzo and E. J. Mele, Phys. Rev. B 29, 1685 
(1984). 

4 N. H. Shon and T. Ando, J. Phys. Soc. Jpn. 67, 2421 
(1998). 

5 Y. Zheng and T. Ando, Phys. Rev. B 65, 245420 (2002). 

6 V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 
146801 (2005). 

7 N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 
Rev. B 73, 125411 (2006). 

8 T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005). 

9 E. McCann and V. I. Falko, Phys. Rev. Lett. 96, 086805 
(2006). 

10 M. Koshino and T. Ando, Phys. Rev. B 73, 245403 (2006). 

11 M. I. Katsnelson, Euro. Phys. J. B 52, 151 (2006). 

12 J. Nilsson, A. H. Castro Neto, N. M. R. Peres, and F. 



Guinea, Phys. Rev. B 73, 214418 (2006). 

F. Guinea, A. H. Castro Neto, and N. M. R. Peres, Phys. 

Rev. B 73, 245426 (2006). 

E. McCann, Phys. Rev. B 74, 161403 (2006). 

E. V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. 

Peres, J. M. B. Lopes dos Santos, J. Nilsson, F. Guinea, 

A. K. Geim, and A. H. Castro Neto, Phys. Rev. Lett. 99, 

216802 (2007). 

M. Koshino, New J. Phys. 11, 095010 (2009). 
K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 
Science 306, 666 (2004). 

K. S. Novoselov, E. McCann, S. V. Morozov, V. I. Falko, 
M. I. Katsnelson, U. Zeitler, D. Jiang, F. Schedin, and A. 
K. Geim, Nature Phys. 2, 177 (2006). 
C. Berger, Z. Song, T. Li, X. Li, A. Y. Ogbazghi, R. Feng, 
Z. Dai, A. N. Marchenkov, E. H. Conrad, P. N. First, and 
W. A. de Heer, J. Phys. Chem. B 108, 19912 (2004). 



12 



20 T. Ohta, A. Bostwick, T. Seyller, K. Horn, and E. Roten- 
berg, Science 313, 951 (2006) 

21 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 
I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. 
Firsov, Nature 438, 197 (2005). 

22 Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
438, 201 (2005). 

23 M. Koshino and T. Ando, Physica E 40, 1014 (2008). 

24 M. Koshino and E. McCann, Phys. Rev. B 81, 115315 
(2010). 

25 M. Fujita, K. Wakabayashi, K. Nakada. and K. Kusakabe, 
J. Phys. Soc. Jpn. 65, 1920 (1996). 

26 K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dressel- 
haus, Phys. Rev. B 54, 17954 (1996). 

27 K. Wakabayashi, Phys. Rev. B 64, 125428 (2001). 

28 E. McCann and V. I. Falko, J. Phys.: Condens. Matter 16, 
2371 (2004). 

29 L. Brey and H. A. Fertig, Phys. Rev. B 73, 235411 (2006). 

30 N. M. R. Peres, A. H. Castro Neto, and F. Guinea, Phys. 
Rev. B 73, 241403 (2006). 

31 K. Wakabayashi, J. Phys. Soc. Jpn. 71, 2500 (2002). 

32 K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, Phys. 
Rev. B 59, 8271 (1999). 

33 K. Wakabayashi and M. Sigrist, Phys. Rev. Lett. 84, 3390 
(2000). 

34 Y.-W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 
97, 216803 (2006). 

35 Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature 444, 
347 (2006). 

36 B. Obradovic, R. Kotlyar, F. Heinz, P. Matagne, T. Rak- 
shit, M. D. Giles, M. A. Stettler, and D. E. Nikonov, Appl. 
Phys. Lett. 88, 142102 (2006). 

37 L. Yang, C.-H. Park, Y.-W. Son, M. L. Cohen, and S. G. 
Louie, Phys. Rev. Lett. 99, 186801 (2007). 

38 L. Yang, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 
101, 186401 (2008). 

39 T. C. Li and S.-P. Lu, Phys. Rev. B 77, 085408 (2008). 

40 H. Raza and E. C. Kan, Phys. Rev. B 77, 245434 (2008). 

41 V. Ryzhii, M. Ryzhii, A. Satou, and T. Otsuji, J. Appl. 
Phys. 103, 094510 (2008). 



42 T. Wassmann, A. P. Seitsonen, A. M. Saitta, M. Lazzeri, 
and F. Mauri, Phys. Rev. Lett. 101, 096402 (2008). 

43 A. R. Akhmerov and C. W. J. Beenakker, Phys. Rev. B 
77, 085423 (2008). 

44 V. H. Nguyen, V. N. Do, A. Bourne, V. L. Nguyen, and P. 
Dollfus, J. Phys.: Conf. Ser. 193, 012100 (2009). 

45 D. Gunlycke and C. T. White, Phys. Rev. B 81, 075434 
(2010). 

46 B. Sahu, H. Min, A. H. MacDonald, and S. K. Banerjee, 
Phys. Rev. B 78, 045404 (2008). 

47 E. V. Castro, N. M. R. Peres, J. M. B. Lopes dos Santos, 
A. H. Castro Neto, and F. Guinea, Phys. Rev. Lett. 100, 
026802 (2008). 

48 D. A. Abanin, P. A. Lee, and L. S. Levitov, Phys. Rev. 
Lett. 96, 176803 (2006). 

49 D. A. Abanin, K. S. Novoselov, U. Zeitler, P. A. Lee, A. 
K. Geim, and L. S. Levitov, Phys. Rev. Lett. 98, 196806 

(2007) . 

50 D. A. Abanin, P. A. Lee, and L. S. Levitov, Solid State 
Commun. 143, 77 (2007). 

51 J. Nilsson, A. H. Castro Neto, F. Guinea, and N. M. R. 
Peres, Phys. Rev. B 76, 165416 (2007). 

52 J. W. Gonzalez, H. Santos, M. Pacheco, L. Chico, and L. 
Brey, Phys. Rev. B 81, 195406 (2010). 

53 T. Nakanishi, M. Koshino, and T. Ando, 
arXiv:1008.4450vl, to be published in Phys. Rev. 
B. 

54 A. Misu, E. Mendez, and M. S. Dresselhaus, J. Phys. Soc. 
Jpn. 47, 199 (1979). 

55 T. Nakanishi and T. Ando, J. Phys. Soc. Jpn. 77, 024703 

(2008) . 

56 Y. Kobayashi, K. Fukui, T. Enoki, K. Kusakabe, and Y. 
Kaburagi, Phys. Rev. B 71, 193406 (2005). 

57 Y. Niimi, T. Matsui, H. Kambara, K. Tagami, M. Tsukada, 
and H. Fukuyama, Phys. Rev. B 73, 085421 (2006). 

58 T. Matsui, H. Kambara, Y. Niimi, K. Tagami, M. Tsukada, 
and H. Fukuyama, Phys. Rev. Lett. 94, 226403 (2005). 

59 Y. Niimi, H. Kambara, T. Matsui, D. Yoshioka, and H. 
Fukuyama, Phys. Rev. Lett. 97, 236804 (2006). 



